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PgQced®, installation <st prograigm© d ' ogdinatQur pour 
es timer 1 ' ef fectif initial d'una population d'acides 
nucleiquaa, notastment par PCR 



La presente invention concerne 1' estimation de I'effectif 
initial d * une population d'interet dans un echantillon 
soumis a succession de reactions d ' amplification . 



La presente invention trouve une application avantageuse 
mais aucunement limitative a la determination d'une 
quantite initiale d'acides nucleiques dans un echantillon 
soumis a une reaction de E^CR (pour "Polymerase Chain 
Reaction") , en temps reel. Une technique de ce type, dite 
"de quantification par PCR" , est utilises notamment pour 
evaluer le nombre de copies d' agents pathogenes (par 
exemple du Virus de 1' Immunodef icience Humaine VIH) dans 
un prelevement de fluides corporels d'un patient, 
typiquement dans le cadre d'un suivi medical. 

On se refere a la figure 1 pour decrire brievement 
1' allure schematique d'une courbe d ' amplification par PCR 
en temps reel, avec, en abscisses, des indices de cycles 
de PCR et, en ordonnees, dans 1' exemple represents, des 
quantites de fluorescence emises (en unite arbitraire) et 
mesurees a chaque cycle de PCR. On indique qu'a chaque 
cycle de PCR, 1' echantillon subit des variations de 
temperature permettant a I'ADN polymerase de proceder a 
1' amplification des acides nucleiques et les produits PCR 
correspondants sont detectes par des molecules 
f luorescentes . On obtient, en reportant la quantite de 
fluorescence mesuree Fn en fonction du nombre n de cycle 
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de PGR, une variation du type represente sur la figure 1 
et comprenant au moins : 

- une premiere partie BDF ou les mesures de fluorescence 
se confondent sensiblement avec un bruit de fond de 

5 I'appareil de mesure de fluorescence, 

- une seconde partie EXP ou les quantites de fluorescence 
mesurees croissenL de fagon sensiblement exponentielle, 

- une troisieme partie LIN ou la croissance des quantites 
de fluorescence mesurees est sensiblement attenuee et 

10 se comporte globalement de fagon sensiblement lineaire, 

et 

- une quatrieme partie PLA ou les mesures de fluorescence 
atteignent une phase de plateau. 

15 On retiendra en fait que pour les premiers cycles de PGR 
(premiere et seconde parties), la population d'interet 
croit sensiblement de fagon exponentielle, tandis que pour 
les cycles qui suivent (troisieme et quatrieme parties), 
d'autres phenomenes entrent en competition avec la 

20 croissance de la population d'interet, cette croissance 
etant alors amortie jusqu'a la phase de plateau PLA. 

On connait notamment par le document : 

"Mathematics of quantitative kinetic PCR and the 
25 application of standard curves", de R.G.Rutledge et 
C.Cote, dans Nucleic Acids Research , 2003, vol.31, N'^IG, 
une methode pour estimer la quantite initiale, inconnue, 
d*acides nucleiques dans un echantillon d'interet, par 
PCR. Elle consists a utiliser plusieurs echantiilons de 
30 quantites initiales connues en acides nucleiques, dits 
"standards" , pour determiner, par interpolation, la 
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quantite initiale d'acides nucleiques presents dans 
1' echantillon d'interet. 

De fagon generale, plus la quantite initiale en acides 
nucleiques dans un echantillon est elevee, et plus tot la 
quantite de produit PGR obtenue et, par consequent, la 
fluorescence emise deviendront detectables. En se referant 
a la figure 2 relative a I'art anterieur, on comprendra 
que la population initiale dans le standard Stl est plus 
elevee que celle du standard St2, laquelle est plus elevee 
encore que celle du standard St3, etc, car le cycle Ctl 
pour le standard Stl intervient avant le cycle Ct2 
correspondant pour le standard St2, lequel intervient 
encore avant le cycle Ct3 du standard St3, etc. 

Ainsi, ce cycle Ct, correspondant au cycle auquel la 
mesure de fluorescence atteint un seuil de fluorescence 
THR (comme represente sur la figure 2), prefixe 
arbitrairement { typiquement au-dessus du bruit de fond), 
se presente comme un parametre representatif de I'effectif 
initial No d'une population d'acides nucleiques soumise 
aux cycles de PGR. Cette observation a ete mise a profit 
dans I'art anterieur cite ci-avant pour etablir une 
dependance telle que representee sur la figure 3 entre les 
numeros de cycle Ctl, Gt2, Gt3, Ct4, pour plusieurs 
standards d'effectifs initiaux connus, et leur effectif 
initial No', Nq% No^ Nq' . En reportant ainsi les cycles 
Gtl, Gt2, Gt3, Ct4, etc, en ordonnees et le logarithme des 
effectifs initiaux NqS No^ No^ No' en abscisses, on 
obtient une pente de regression REG. Sur cette pente REG, 
on reporte (fleche Fl en traits pointilles) le cycle Gtint 
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detecte pour 1 ' echantillon d'interet. Par interpolation 
sur la pente de regression REG (fleche F2 en traits 
pointilles), on determine ensuite I'effecrif initial Nq^"^*^ 
de la population dans 1' echantillon d'interet. 

5 

Cette methode, meme si elle s'est largement repandue, 
presente neanmoins quelques inconvenieiits . 

Tout d'abord, elle necessite 1 ' utilisation de plusieurs 
10 echantillons standard d*effectifs initiaux respectifs 
connus . 

Ensuite, la methode depend du juge de 1 ' ut ilisateur car la 
valeur du seuil de fluorescence qui est choisie par 

15 1' utilisateur influence directement la valeur des cycles 
Ct des courbes d' amplification et, par consequent, la 
valeur estimee de I'effectif initial dans 1 ' echantillon 
d'interet. La valeur du seuil a egalement un impact sur la 
precision du resultat car la precision est generalement 

20 raeilleure si le seuil est choisi dans la phase EXP de 
croissance exponentielle des courbes d' amplification , 
Toutefois, en pratique, il est difficile pour 
1' utilisateur de savoir si le niveau seuil de fluorescence 
fixe THR correspond bien a la phase exponentielle des 

25 courbes, et ce, pour tous les echantillons (les 
echantillons standard et 1 * echantillon d'interet). 

Enfin, la methode sous-entend, sans aucune verification, 
que la population a le meme rendement d' amplification dans 
30 1 ' echantillon d'interet et dans tous les echantillons 
standard. Ainsi, si typiquement 1' echantillon d'interet 
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contient des inhibiteurs de PGR, son resultat sera 
f aussement abaisse . 

On retiendra done que la technique de I'art anterieur 
depend du seuil de fluorescence THR defini par 
1 ' utilisateur . La valeur choisie influence les valeurs des 
cycles Ct et par consequent la determination de la 
quantite initials dans 1 ' echantillon d'interet. C'est 
I'une des raisons pour lesquelles de nombreux travaux ont 
cherche recemment a automatiser et a fiabiliser la 
detection du cycle Ct . 

La presente invention vient ameliorer la situation en 
proposant une approche completement differente. 

Elle vise tout d'abord un precede, ce procede etant mis en 
oeuvre par des moyens inf ormatiques pour quantifier, de 
fagon absolue et/ou relative, une population initiale 
d'acides nucleiques dans un echantillon d'interet. Get 
echantillon est soumis a une succession d • applications 
d'une reaction d ' amplification de la population d'interet. 
De fagon tres generale, cette amplification peut etre 
menee par cycles successifs de PGR, mais aussi par toute 
autre technique d ^ amplification . On retiendra surtout que 
cette amplification doit simplement etre definie par un 
rendement de reaction, comme on le verra plus loin. Au 
cours de ces amplifications successives, on releve des 
mesures experimentales representatives d'un effectif 
courant de la population au moins dans 1 ' echantillon 
d'interet. On comprendra que 1 ' on peut relever une ou 
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plusieurs mesures apres ou pendant chaque reaction 
d' amplification, sans parte de generalite. 

Salon una definition actuellement preferee de 1* invention, 
le procede au sens de 1' invention comporte les etapes 
suivantes : 

a) prevoir un modele de rendement de la reaction 
d ' amplification en fonction de la succession 
d' amplifications, ce modele comportant : 

o una phase sensiblement constante pour une premiere 

partie des applications de la reaction, et 
o une phase non constante pour une seconde partie 
des applications de la reaction, 
les premiere et seconde parties etant separees par une 
region de basculement du rendement entre les phases 
constante et non constante, incluant un indice 
d ' application correspondant sensiblement au 

basculement , 

b) exprimer, en utilisant le modele de rendement, une 
relation faisant intiervenir au moins 1 ' indice de 
basculement et un parametre representant I'effectif 
initial dans 1 ' echant illon d'interet, 

c) par comparaison avec les mesures experimentales , 
determiner au moins 1 ' indice de basculement, et, dans 
une etape d) ulterieure ou conjointe, en deduire 
I'effectif initial dans 1 ' echantillon d'interet. 

D'autres avantages et caracteristiques de 1' invention 
apparaitront a la lecture de la description detaillee du 
mode de realisation donne ci-apres a titre d' example, en 
s*appuyant sur les graphiques annexes et sur lesquels : 
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outre la figure 1, relative a I'art anterieur, 
representant une variation de la quantite de 
fluorescence mesuree en fonction d'un nombre de cycles 
de PGR, decrite ci-avant ; 

outre la figure 2, relative a I'art anterieur, 
representant les quantites de fluorescence emises et 
croissantes en fonction du nombre de cycles de PGR, et 
decrite ci-avant ; 

outre la figure 3 representant une methode 
d ' interpolation pour determiner la quantite initiale de 
la population d'interet dans 1 * echantillon d'interet, 
selon une methode connue de I'art anterieur, et decrite 
ci-avant ; 

la figure 4A represente une variation schematique des 
mesures exper imentales precitees en fonction de la 
succession d ' amplifications que 1 • on applique a 
1 ' echantillon d * interet , 

la figure 4B represente une variation schematique du 
rendement de la reaction d ' amplification, obtenue a 
partir des mesures experimentales , en fonction de la 
succession d ' amplifications appliquees a 1 ' echantillon 
d' interet , 

la figure 5 represente une regression entre des indices 

de basculement du rendement entre la phase constante et 
la phase non constante et le logarithme des effectifs 
initiaux pour des echantillons standard et pour 
1' echantillon d'interet, pour la mise en oeuvre d'un 
premier mode de realisation, 

la figure 6A represente une variation typique de la 
quantite de fluorescence mesuree et ajustee en tenant 
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compte d'un bruit de fond propre aux mesures, en 
fonction du nombre de cycles n de la reaction PCR, 

- la figure 6B represente une variation de I'efficacite 
de la PGR de la figure 6A, en fonction du nombre de 

5 cycles n, 

- la figure 7 represente une comparaison de la variation 
experimentale de la fluorescence emise selon la figure 
6A et le resultat obtenu par 1 ' application du modele de 
la fluorescence emise obtenu par integration d'un 

10 modele d ' ef f icacite, selon un second mode de 

realisation, 

- la figure 8 represente une comparaison de la variation 
de I'efficacite de la figure 6B avec 1 ' application d'un 
modele d ' ef f icacite, deduit de la modelisation de la 

15 fluorescence emise de la figure 7, 

- la figure 9 represente 1 ' organigramme schematisant les 
principales etapes du precede selon le mode - de 
realisation propose pour la presente invention, et 

la figure 10 represente une installation pour 
20 quantifier la population initiale d'un echantillon 

d ' interet . 

On s'appuie sur les figures 4A et 4B pour decrire 
brievement quelques principes de 1' invention illustrant 
25 les caracteristiques du procede ci-avant. 

Tout d'abord on indique que la figure 4A represente une 
succession de mesures experimentales Fn representatives de 
I'effectif courant d'une population d'interet qui subit 
30 progressivement une succession de reactions 

d ' amplification, chaque reaction etant indexee par un 
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indice n. Cette succession de reactions correspond, dans 
1' example decrit ici a titre non limitatif, a une 
succession de cycles de PGR. Les mesures exp§rimentales Fn 
correspondent dans cet exemple a titre non limitatif a des 
mesures de quantites de fluorescence a chaque cycle de 
PGR. Ainsi, dans le cadre d ' une methode de quantification 
qui combine la reaction PGR et la fluorescence emise par 
I'echantillon d'interet, on introduit des reactifs 
fluorescents dans I'echantillon, de sorte que la 
fluorescence emise, au cours d'un cycle de PGR, est 
proportionnelle a I'effectif d'acides nucleiques dans cet 
echantillon. Neanmoins, cette realisation admet de 
nombreuses variantes. Deja, il peut etre prefere de mener 
plusieurs mesures ou aucune mesure pour certains cycles de 
PGR. Par ailleurs, plus generalement , la methode de mesure 
peut utiliser d'autres techniques que la fluorescence, 
meme si cette derniere est souvent celle employee dans la 
quantification par PGR. Enfin, on indique que d'autres 
techniques d ' ampl if ication peuvent etre prevues dans le 
cadre de la presente invention, des lors qu'il est 
possible de suivre la variation du rendement de la 
reaction correspondant a cette amplification. Neanmoins, 
comme 1' exemple decrit ci-apres traite pref erentiellement 
de cycles de PGR, on parlera d' " efficacite de la PCR" , 
notee En pour chaque cycle de PGR d' indice n, pour 
designer le rendement de la reaction d ' amplification . 

Gomme indique ci-avant en reference a la figure 1, on 
retrouve sur la figure 4A principalement deux regions ou : 
- pendant les premiers cycles de PGR (partie EXP) , la 
population croit sensiblement de fagon exponentielle , 



1 er depot 



10 

- tandis que pour les cycles qui suivent (partie LIN et 
PLA) , d'autres phenomenes entrent en competition avec 
la croissance de la population d'interet, cette 
croissance etant alors amortie. 

5 

On formule alors les hypotheses suivantes : 

- le rendement de la reaction En est relativement constant 
durant les premiers cycles sur la partie EXP, et 

- a partir d ' un certain nombre de cycles effectues, le 
10 rendement En de la reaction decroit sur la partie LIN et 

PLA. 

Cette decroissance du rendement peut avoir diverses 
explications, notamment une degradation et/ou une carence 
des reactifs de PGR (DNA polymerase, dNTPs, amorces, ou 
15 autres) et/ou une inhibition par les produits formes eux- 
memes . 

On fait ainsi I'hypothese ici que le rendement est dVabord 
constant et decroit ensuite* Toutefois, on comprendra que 
20 1' invention s* applique de fagon plus generale dans le 
cadre d'un rendement : 

- d*abord constant, ce qui correspond a une situation 
normale de croissance par amplification, 

- puis non constant (decroissant ou croissant), ce qui 
25 correspond a une situation sensiblement anormale. 



On a constate que dans le cadre de reactions 
d ' amplification de quantite d'acides nucleiques, le 
rendement bascule souvent d'une phase constante a une 
30 phase non constante. Cette observation a ete mise a profit 
au sens de 1* invention pour en deduire la quantite 
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initiale d'acides nucleiques, comme on le verra en detail 
ci-apres. Auparavant, on indique simplement que le 
rendement peut encore basculer d'une phase non constante 
pour les premiers cycles a une phase constante ensuite. La 
presente invention s» applique tout aussi bien a ce cas . De 
maniere generale, on retiendra done que 1 ' on detecte, au 
seiis de I'invenLiou, un basculement du rendemenL, d'une 
phase constante a une phase non constante. 

On rappelle que le but recherche est de donner I'effectif 
initial de la population soumise a amplification. Deja, en 
se referant a la figure 4A, on comprendra que la mesure Fo 
representative de cet effectif initial, pratiquement 
confondue avec le bruit de fond des mesures BDF ne permet 
pas directement d'acceder, a elle seule, a I'effectif 
initial. Dans I'art anterieur, on a recherche a quantifier 
cet effectif initial en exploitant la phase exponentielle, 
done typiquement a la sortie du bruit de fond. On y 
determine alors le cycle seuil Ct (correspondant au point 
AA pour "art anterieur''} de la figure 4A. Or, comme on I'a 
vu ci-avant, les mesures sont souvent bruitees dans cette 
region et il est difficile de determiner avec precision un 
cycle seuil Ct de sortie du bruit de fond. 

Dans une toute autre approche, la presente invention 
expioite plutot la quasi-totalite des points de la courbe 
d ■ amplification pour determiner precisement une region BAS 
ou le rendement bascule entre une phase constante et une 
phase non constante, typiquement ici entre la phase 
exponentielle EXP et la phase lineaire LIN. Deja, on 
constatera que les mesures sont logiquement moins bruitees 
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dans cette region BAS que dans la region de sortie du 
bruit de fond puisque la region BAS intervient pour des 
cycles ulterieurs. De plus, notamment en raison des 
proprietes mathematiques liees au rendement, on verra plus 
5 loin que le nombre de standards qui peuvent etre utilises 
pour la quantification de I'effectif initial dans la 
population d'inLeret est, de fagon tres avantageuse, 
reduit par rapport au nombre de standards utilises dans la 
quantification de 1 ' art anterieur. 

10 

On decrit maintenant brievement la relation qui peut lier 
la region de basculement BAS a I'effectif initial de . la 
population d'interet. On rappelle au prealable que ie 
rendement d-une reaction d ' amplification est donne par : 

15 Nn+l = Nn + En.Nn, OU : 

- Nn est I'effectif de la population d'interet apres une 
amplification d'indice n dans une succession 
d' amplifications, 

- Nn+i est I'effectif de la population d'interet apres une 
20 amplification suivante, d'indice n+1, dans la 

succession d ' amplifications precitee, et 

- En est le rendement de la reaction d ' amplification 
d'indice n dans la succession d ' amplifications 
precitee . 

25 En reformulant cette relation, par recurrence, on 

obtient : 

Nn+l = (1 + En) (1 + En-l) d + En-2) - d + Eq) Nq, 

ou No est I'effectif initial de la population d'interet. 
Tant que le rendement En est constant, on comprendra que 
30 cette derniere relation peut s'ecrire simplement : 
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Nn-fi = (1 + Eo)''^^ . No (ou I'indice n + 1 n'a pas 
encore atteint la region de basculement BAS) . 
En effet, le rendement etant alors constant pour ces 
premiers cycles, on a : 

5 En = En-l ~ En-2 = ••• — Eo/ 

ou Eo est la valeur de la phase constante du rendement. 
Neanmoins, lorsque 1 * indice n + 1 depasse la region de 
basculement BAS, la relation devient : 

Nn+i = No X (1 + Eo)*^EEP X fonction (CEEPr n-hl) , ou : 
10 - (Ceep - 1) est le dernier indice de la reaction 
d ' amplification au cours duquel le rendement est encore 
constant (on comprendra done que 1' indice Ceep 
represente lui-meme 1 ' indice de basculement, proprement 
dit, entre la phase de exponentielle a la phase 
15 lineaire) , 

- le terme fonction (Ceep^ n+1) est une fonction 
particuliere qui caracterise la phase non constante du 
rendement et qui depend au moins de 1' indice de 
basculement Ceep et de 1' indice courant de 
20 1 ' amplification n+1 . 



On voit ainsi comment peuvent etre lies 1 ' indice de 
basculement Ceep et I'effectif initial No de la population 
d'interet. On comprendra qu ' a ce stade, on a deja realise 
25 les etapes a) et b) du precede tel que defini ci-avant. 



Deja, un premier mode de realisation consiste a determiner 
experimentalement 1 ' indice de basculement Ceep et a le 
correler a I'effectif initial, par regression, en 
30 utilisant plusieurs echantillons standard qui subissent le 
meme traitement d ' amplification que 1 ' echantillon 
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d'interet. Dans ce cas, on comprendra que les etapes b) et 
c) du procede tel que defini ci-avant sont simplement 
inversees puisque 1 ' on determine experimentalement , 
d'abord, I'indice de basculement Ceep (etape c) ) et que 
5 I'on determine ensuite la relation entre 1 ' indice Ceep et 
I'effectif initial No (etape b) ) pour obtenir finalement 
I'effectif initial No (etaped)). 

Avant de decrire en detail toutes ces etapes au sens de ce 
10 premier mode de realisation^ on decrit un procede de 
determination de 1 ' indice Ceep/ a partir des mesures 
experimentales . En particulier, on comprendra que cette 
maniere de determiner experimentalement 1 ' indice Ceep peut 
etre appliquee a un autre mode de realisation, different 
15 du premier mode de realisation precite. 

En reprenant la relation entre 1 ' ef f icacire En a un. cycle 
donne n et I'effectif courant dans 1 ' echantillon d'interet 
au meme cycle Nn et a un cycle ulterieur Nn+i, on peut 
20 exprimer I'efficacite de 1 ' amplification comme suit : 
En = Nn..i/N„ - 1 

Dans certains cas notamment ou il n*est pas necessaire de 
tenir compte du bruit de fond BDF dans les mesures, on 

25 peut en premiere approximation considerer que les mesures 
sont sensiblement proportionnelles a I'effectif courant 
dans 1 ' echantillon d'interet. Toutefois, en pratique, on 
tiendra le plus souvent compte des derives des mesures et, 
dans une realisation qui sera decrite maintenant, on 

30 determine des mesures experimentales corrigees, notees 
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F'nf a partir des mesures directes Fn telles que 
representees sur la figure 4A. 

Preferentiellement , on applique une etape prealable de 
traitement des mesures experimentales Fn consistant a 
soustraire le bruit de fond BDF et a introduire ensuite 
une conperisation pour tenir compte d'une mesure non nulie 
e representative de l^effectif initial. Dans I'exemple 
represents sur la figure 4A, la variation du bruit de fond 
BDF en fonction de I'indice n peut etre representee par 
une fonction lineaire car les essais ont montre qu ' un 
modele lineaire est satisfaisant pour des mesures de 
fluorescence en PGR. Neanmoins, un modele de variation 
exponentielle pourra etre prefere dans certains cas. On 
applique en tout etat de cause le modele qui respecte le 
mieux la variation du bruit de fond BDF donnee typiquement 
par les premiers points de mesure. Ensuite, on soustrait 
le modele BDF de variation de bruit de fond choisi a 
toutes les valeurs de mesures experimentales Fn. En 
appliquant cette etape, on comprendra que la mesure de 
fluorescence theorique Fo est ramenee a une valeur de 
mesure nulle, ce qui correspond a un effectif initial No 
nul, ce qui n'a aucune realite physique. Par consequent, 
on applique avantageusement une compensation de cette 
correction comme suit : 

F*n = - BDF + ou : 

- le terme F'n correspond a une mesure corrigee pour un 
indice courant n, 

- le terme Fn correspond a la mesure experimentale brute a 
cet indice courant n. 
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- le terme BDF correspond a la valeur du bruit de fond 
modelise pour I'indice courant n, et 

- € est le terme de compensation correspondant , suppose 
constant dans 1 ' example decrit et representant 

5 directement en fait I'effectif initial No. 

Ces etapes de correction par rapport au bruit de fond, 
meme si elles sont tres avantageuses pour la determination 
de I'indice de basculement C EEPf peuvent etre appliquees 

10 par ailleurs a toute determination et quantification de 
I'effectif initial No des lors qu ' un bruit de fond est 
susceptible de fausser la mesure de cet effectif No- A ce 
titre, ces etapes peuvent faire I'objet d'une protection 
separee, le cas echeant- 

15 

Les . mesures corrigees F*n ainsi obtenues sont 
avantageusement proportionnelles aux effectifs courants Nn 
dans 1 ' echant illon d'interet, de sorte que le rendement En 
peut maintenant s'exprimer directement en fonction des 
20 valeurs de mesure (corrigees comme indique ci-avant) , par 
la relation : 

En = E'n+i/F*n ~ 1 

Ainsi, a partir des mesures experimentales Fn de la figure 
25 4A, on obtient les mesures corrigees F'n a partir 
desquelles on determine ensuite la variation de 
I'efficacite En en fonction de la succession d'indices n, 
telle que representee sur la figure 4B. 

30 En bref, on exprime les mesures experimentales sous la 
forme d*une variation experimentale du rendement En telle 
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que representee sur la figure 4B, en fonction de la 
succession d ' applications n. On obtient alors une 
variation exper imentale du rendement, comportant : 

- une premiere region sensiblement bruitee pour des 
indices d ' applications n faibles {sensiblement avant le 
cycle CG dans 1 ' exemple de la figure 4B) , 

- ' suivie d'une seconde region moins bruitee poxir des 

indices d ' applications plus eleves (au moins apres la 
region de basculement BAS) . 

Au moins dans le cas le plus courant en amplification par 
PGR et en mesure par fluorescence, la phase non constante 
du rendement est decroissante et correspond a cette 
seconde region peu bruitee (comme represents sur la figure 
4B) - Notamment pour eliminer des points de mesure qui 
risqueraient de fausser les resultats lorsque 1 ' on choisit 
un modele a appliquer a la variation du rendement : 

- on estime une valeur grossiere de la phase constante du 
rendement Eq, et 

- notamment pour la recherche de I'indice de basculement 
Geep/ on ignore une partie au moins des mesures dans la 
seconde region moins bruitee et dont le rendement 
estime est inferieur a une valeur seuil, telle que, par 
exemple, une fraction de la phase constante Eq. 

Ces points NEG (figure 4B) qui sent elimines sont 
typiquement ceux cor respondant a des indices 
d ' ampl:^f ication n tres eleves et qui pourraient ne plus 
respecter le modele qui a ete choisi pour I'efficacite 
sensiblement autour de la region de basculement BAS. Pour 
les eliminer, on evaluera par exemple une moyenne de la 
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phase constante du rendement Eq typiquement pour les 
premiers indices n. Ensuite, on choisira une valeur seuil 
correspondant a une fraction, par exemple 10%, de la 
moyenne trouvee pour la phase constante Eo- Ensuite, en 
5 partant des indices n les plus eleves, on elimine tous les 
points de mesure NEG dont la mesure de rendement est 
inferieure ou eyale a cette valeur seuil. Cette etape, 
tres avantageuse pour la detection de 1 ' indice Ceep, peut 
neanmoins etre appliquee a toute determination bases sur 
10 le rendement En et, le cas echeant, peut faire I'objet 
d'une protection separee. 

Dans le cas d'une phase non constante du rendement qui est 
decroissante et qui suit la phase constante, comme 

15 represents sur la figure 4B, on repere la region de 
basculement BAS, dans le sens des indices n decroissants a 
partir de la seconde region moins bruitee, en detectant un 
indice grossier CG dont le rendement depasse une valeur 
predeterminee . Ainsi, en se referant a la figure 43, on va 

20 dans le sens des indices n decroissants en remontant vers 
la region de basculement BAS et on evalue le rendement 
associe a cet indice. Pour le premier point de mesure dont 
le rendement depasse signif icativement la valeur 
predeterminee precitee, on considere que 1' indice grossier 

25 precite CG est detecte et correspond a 1 ' indice de ce 
point de mesure. 



Comme on le verra dans un 
on peut modeliser, pour 
30 variation de son rendement 
ef f ectivement a 1* indice 



mode de realisation ulterieur, 
chaque point de mesure, la 
comme si ce point correspondait 
de basculement Ceep- Dans cette 
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realisation, on estime alors sa phase constante Eq de 
rendement et si cette valeur estimee depasse la valeur 
predeterminee precitee, on considere ce point conune 
correspondant a 1 ' indice grossier CG. 

De fagon generale, un rendement maximum a une valeur de 1 
de sorte que 1 * on peut choisir la valeur predeterminee 
precitee comme etant egale a 1. Neanmoins, ce choix admet 
des variantes et par exemple, il peut etre prevu de fixer 
cette valeur predeterminee comme correspondant a la 
moyenne du rendement Eo evaluee sur les premiers cycles de 
reaction . 

Ensuite, on affine 1' estimation de la valeur, qui peut 
avantageusement etre f ractionnaire, de I'indice 
d' application Cee? dans la region de basculement, dans le 
sens des indices d ' applications croissants a partir de 
I'indice grossier CG, en detectant un indice d ' application 
d' amplification dont le rendement vaut environ la valeur 
predeterminee precitee. Ainsi, en se referant a nouveau a 
la figure 4B, pour af finer la recherche de I'indice de 
basculement Ckkp apres la determination de I'indice 
grossier CG, on redescend a partir de I'indice grossier CG 
dans le sens des indices n croissants, par pas plus petit 
qu'un indice entier et on determine, par exemple par 
interpolation, I'abscisse pour laquelle on franchit la 
valeur predeterminee. Typiquement, tant que la valeur de 
phase constante reste superieure a 1, on poursuit la 
recherche dans le sens des n croissants et 1 ' on determine 
I'indice (Ceep - D qui precede le basculement des que la 
valeur constante Eo est egale ou est tres proche de 1. 
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C'est la raison pour laquelle on choisira un pas de 
recherche cor respondant a une fraction d'indice, par 
example 10% d ' un cycle n, 

Dans le premier mode de realisation precite, on prevoit 
une pluralite d ' echantillons standard d'effectifs initiaux 
respectifs connus, auxquels on applique la" succession 
d ' amplifications , sensiblement dans les memes conditions 
que pour 1 ' echantillon d'interet. On determine aussi leurs 
indices respectifs de basculement selon les etapes a) , b) 
et c) decrites precedemment . A 1 * etape d) : 

- on etablit une dependance entre les effectifs initiaux 
des echantillons standard Nq^^ et leurs indices Cs^p^^, 

- et, apres la determination de 1 * indice pour 
1 ' echantillon d*interet Ceep/ on deteirmine l^effectif 
initial de la population d'interet No par interpolation 
sur cette dependance. 

Ainsi, en se referant a la figure 5, on peut etablir une 
dependance entre le cycle de basculement Ceep^/Ceep^ des 
standards et leur concentration initiale Nq^/No^ (en fait 
leur logarithme) , par exemple par regression. En mesurant 
I'indice de basculement Ceep^''^ pour 1 ' echantillon d'interet 
et en reportant sa valeur sur la pente de regression de la 
figure 5, on obtient par interpolation la concentration 
initiale No^"^ dans 1 ' echantillon d'interet. 

Ce premier mode de realisation s ' apparente done 
sensiblement a celui de I'art anterieur decrit en 
reference a la figure 3. Cependant, on gardera a 1' esprit 
que I'indice de basculement Ceep sur lequel se base ce 
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premier mode de realisation ne correspond aucunement au 
cycle Ct de I'art anterieur, 

Dans une approche sensiblement differente de ce premier 
mode de realisation : 

- a I'etape b) , on exprime en fait, en utilisant le 
modele de rendement, une varialion paraiaetree en 
fonction de la succession d * amplifications , cette 
variation mettant en jeu au moins un parametre 
representant I'indice de basculement Ceep/ ^t 

- a I'etape c) , par comparaison avec les mesures 
experimentales, on determine au moins ce parametre 
representant I'indice de basculement Ceep- 

Dans un second mode de realisation, cette variation 
parametree est representative de I'effectif courant Nn 
dans 1 • echantillon d'interet. 

Typiquement, cette variation parametree peut etre tiree de 
1' expression du type donne ci-avant : 

= No X (1 + Eo)^EEP X fonction (Ceep^ 
Ainsi, cette variation met en jeu, outre un parametre 
representant 1' indice de basculement Ceep/ un parametre 
representatif de I'effectif initial dans 1 • echantillon 
d'interet Nq. 

Ensuite, aux etapes c) et d) de ce second mode de 
realisation, on determine sensiblement con j ointement ces 
deux parametres Ceep et No- 
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Auparavant, a I'etape a), il faut determiner un modele 
pour la fonction precitee fonction (CsBPf n-i-1) . 

Le plus souvent, pour une quantification par PGR, on 
5 choisira un modele pour la phase non constante du 
rendement correspondant a une exponentielle decroissante 
comportant un parametre de • decroisyance 13 qui sera plus 
detaille ci-apres. Ce parametre de decroissance p est 
alors determine, a I'etape c) , au moins avec I'indice de 
10 basculement Ceep/ P^^ comparaison avec les mesures 
experimentales . 

Ainsi, dans ce second mode de realisation, une fois le 
modele du rendement En choisi, on 1' applique a 
15 1' expression generale de I'effectif courant Nn donnee par 
la relation ci-avant . On obtient alors un modele de la 
variation de I'effectif courant Nn. 

Toutefois, a moins que les mesures experimentales ne 
20 donnent directement la valeur de I'effectif courant Nn (ce 
qui est rarement le cas en pratique, actuellement ) , il 
convient de modeliser ensuite les mesures experimentales 
Fn, elles-memes, en gardant a 1' esprit la soustraction du 
bruit de fond et la compensation consecutive e decrites 
25 ci-avant. 

Ainsi, dans une realisation actuellement preferee, la 

variation parametree precitee : 

- est representative des mesures experimentales, et 
30 - comporte un parametre correspondant a une valeur de 
mesure Fq representative de I'effectif initial. 
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On determine ensuite la valeur de mesure de I'effectif 
initial Fo par comparaison de cette variation parametree 
Fn avec les mesures experimentales . 

En fait, pour mener cette comparaison, on pourra par 
exemple ajuster les parametres Fq, Eq, Ceep et le parametre 
de decroissance p dans le modele des mesures Fn par 
correlations statistiques (typiquement la methode des 
moindres carres) qui s'appuient sur les mesures 
experimentales brutes. Un exemple de realisation est 
detaille ci-apres . 

On obtient d' abord une variation de la quantite de 
fluorescence mesuree et ajustee, en fonction du nombre de 
cycles de PGR appliques, telle que representee par exemple 
sur la figure 6A. Cette figure montre la courbe 
d' amplification d'un echantillon d' inter§t contenant des 
acides nucleiques, ici un fragment d'ADN de quantite 
initiale de 100000 copies, marque par 1 ' intercalant 
SYBRGREEN lors de la reaction de PGR que 1 ' on effectue sur 
I'appareil I-Cycler IQ® de la societe BIO-RAD®. 

Dans 1' exemple decrit, on comprendra que la reaction 
d' amplification est une reaction de PGR en temps reel. Les 
mesures experimentales representent des quantites de 
fluorescence emises. 

On note ci-apres Fn la fluorescence au cycle n 
prealablement ajustee pour le bruit de fond, comme indique 
ci-avant. On note aussi Fo la fluorescence initiale 
theorique avant le premier cycle. On note E„ I'efficacite 
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de la PGR au cycle n. On designe par N le nombre total de 
cycles effectues durant la reaction de PGR. 

Par hypothese, la fluorescence mesuree a chaque cycle n de 
5 la reaction de PGR est definie par : 

Fn.i ^ Fn(l+E„) pour tout n€{0,l,2,...,N-l} (1) 
avec 0 ^ En ^ 1 . 

On calcule I'efficacite de la reaction a chaque cycle n, 
JO comme suit : 

E * ^iitl . 1 pourtoutne{l,2,...,N-l} (2) 

On remarquera que 1' equation (1) est supposee vraie pour 
n==0 . Toutefois, par definition, la fluorescence initiale 
15 Fo est inconnue- L'efficacite au premier cycle Eo ne peut 
done pas etre calculee directement par la formule (2) . 

On a represente sur la figure 6B I'efficacite de la 

reaction PGR approximee par la formule (2) et a partir de 
20 la variation de la fluorescence ajustee de la figure 6A^ 
en fonction du nombre de cycles n. 

On formule pref eren tiellement les hypotheses suivantes : 

- l*efficacite de la reaction est relativement constante 
25 durant les premiers cycles, et 

- a partir d*un certain nombre de cycles effectues, 
I'efficacite de la reaction decroit. 

La figure 6B vient conforter la seconde hypothese, puisque 
30 1 ' on peut constater la decroissance de I'efficacite a 
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partir du cycle n=17 . En revanche, les efficacites des 
cycles 1 a 16 sont fortement bruitees ce qui rend 
difficile la verification graphique de la premiere 
hypothese . 

Toutefois, on affecte pref erentiellement a la variation de 
I'efficacite un niodele du type incluant : 

- une premiere phase, constante, entre le premier cycle 
de PGR et le cycle (Ceep-1) qui precede le cycle de 
basculement note Ceep/ et 

- une seconde phase, decroissante, pour les cycles 
superieurs ou egaux au cycle (Ceep~1) . 

Le cycle (Ceep-1) represente done le dernier cycle, qui 
peut etre f ractionnaire, auquel I'efficacite est encore 
constante . 

On propose alors de modeliser I'efficacite de la reaction 
de la fagon suivante : 



Eo pour 0<n <(C,^^-\) 

(1 \t^^P {-P(n-C,r,.p + 1)) - (3) 

1 + Eo) -I pour (C,>,p-l)<n < (N-1) 



ou Eo et p sont des parametres reels qui sont estimes 
grace a la courbe d ' amplification de la figure 6A, ou 
encore grace a la courbe d'efficacite de la figure 6B, 
comme on le verra ci-apres. 

En variante, on indique qu ' un autre choix, par exemple 
parmi les modeles ci-apres Fl a F3, peut etre prefere, 
notamment en fonction du type d'acide nucleique que 1 ' on 
souhaite quantifier : 
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Fl : 


E„ 


= exp (-p(n-C,,p+l)) - 1 


F2 : 


E„ 


= exp (-n(n-CEEP+l)") - 1 


F3 : 


E„ 


= a - exp (-n(n-CEEp+ir) - 1 



On estime pref erent iellement plusieurs jeux de parametres 
a I'etape c) pour plusieurs candidats de cycle de 
5 basculement Ceep/ et 1 ' on choisit le cycle candidat minimal 
dont les parametres associes maximisent les correlations 
statistiques qui peuvent etre menees a I'etape c) , en tant 
que cycle de basculement Ceep- 



10 Comme on I'a vu precedemment ^ 1* expression (1) ci-avant 
peut encore s'ecrire sous la forme : 

Fn = Fo flO^EJ Pourn € {l,2,...,N} (4) 

Ainsi, en introduisant I'expression de I'efficacite (3) 
15 dans la formule (4), on obtient un nouveau modele a quatre 
parametres (Fo, Eo, p/ Ceep) pour la fluorescence emise 
ajustee ' Fn : 

_ jpoO+Eo)" pour l<n < C,,p 

Fo(l+Eo) ^'-p(P)-« pour Cj^pp < n < N 

20 L'effectif initial dans la population d'interet Nq, 
I'efficacite Eo de la reaction a n = 0, le parametre (3 et le 
cycle de basculement Ceep/ sent evalues repet i t i vement pour 
plusieurs valeurs de cycles dans la region de basculement 
BAS, pour trouver un maximum de correlations statistiques 

25 atteint pour une valeur de cycle minimale egale au cycle 
de basculement Ceep- 
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Dans ce second mode de realisation, on prefere modeliser 
la variation des quantit^s de fluorescence mesurees et 
ajustees, en fonction des cycles, en se basant sur la 
mod^lisation de la variation de 1 ' ef f icacite, et mener 
ensuite les correlations directement sur les quantites de 
fluorescence mesurees et ajustees. 

On remarquera qu ' en ajustant les mesures de fluorescence 
emises pour le bruit de fond, on cree aussi, de fagon 
factice, un ajustement de la fluorescence initiale Fq. 
Ainsi, estimer les parametres du modele d'efficacite ^ 
partir de mesures d'efficacite qui sont d^duites des 
mesures de fluorescence ajustees, est une source d'erreurs 
supplementaires et il peut etre prefere de proceder en 
deux etapes, comme on le verra dans le troisieme mode de 
realisation qui sera decrit plus loin. 

Toutefois, le second mode de realisation propose est plus 
simple et s'adapte bien aux quantifications par PGR 
utilisant des mesures de fluorescence. II se base sur les 
mesures r§elles de fluorescence F'„, qui correspondent aux 
mesures de fluorescence ajustees pour la derive du bruit 
de fond avec une compensation e men§e sur ces mesures. Le 
bruit de fond ayant ete soustrait, on retrouve une 
relation du type : 

ou s est une quant ite qui peut dependre ou non du cycle n. 
Elle est choisie ici constante, pref erentiellement . 
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Dans ce cas, I'efficacite mesuree et "ajustee" aussi, 
notee E'n cycle n est definie par : 



= - 1 = -^i^l-E . 1 pour tout ne{l,2,...,N-l} 



FV 



(7) 



Le modele de la relation (5) ci-avant devient 



F'o(l+E'o)" -E 



pour 1 < n < C 



EEP 



Po(l+E'o)^"'' -£ pour C£Ep< n < N 



- exp (-P'Cn-CHtrp )) 



(8) 



Dans ce cas, les valeurs d'efficacite E'n sont approximees 
experimentalement a partir des mesures, de fagon a pouvoir 

10 fixer un seuil d'efficacite minimal tolere dans la phase 
d'efficacite decroissante . On determine ainsi un cycle 
seuil au-dela duquel les mesures de fluorescence ajustees 
(points NEG figure 4B) ne seront pas exploitees pour mener 
la modelisation , Typiquement, le cycle seuil correspond au 

15 premier cycle de la phase d'efficacite decroissante pour 
lequel I'efficacite est inferieure au seuil d'efficacite 
minimal tolere (par exemple O^l.Eo). 



Plus generalement , la valeur du seuil d'efficacite peut 

20 pref erentiellement etre comprise entre 0 et 0,5, et une 

efficacite de PGR inferieure a cette valeur est 

potentiellement biaisee par des phenomenes d' inhibition 
incontroles - 



25 Dans 1* exemple represents sur la figure 1, la valeur seuil 
d'efficacite a ete fixee a 0,02 (soit 2% de Eq) . Le cycle 
seuil Cs correspondrait au cycle n=36. On a represente sur 
cette figure 7 les mesures de fluorescence emises 
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ajustees. On constate une correlation satisf aisante du 
modele (trait plein) avec les mesures experimentales 
(marquees "o") jusqu'au cycle n=36. La figure 8 montre 
aussi une bonne correlation avec les mesures 
experimentales des predictions de 1 ' ef f icacite , obtenue a 
partir de la figure 7 et de la modelisation de la 
fluorescence mesuree et ajustee. 

Les etapes principales de ce mode de realisation peuvent 
se resu.-ner comme suit, en reference a la f:.gure 9. 

A I'etape de debut 70, on a obtenu les valeurs de 
quantites de fluorescence mesurees et ajustees par rapport 
au bruit de fond en fonction du nombre de cycles n, comme 
represente sur la figure 6A- 

A I'etape 71, on calcule une approximation de I'eff icacite 
de la reaction au cycle n grace a la formule (2) ci-avant, 
pour tout cycle n = 1, 2, (N-1). 

A I'etape 72, on determine le cycle Cs minimal de fagon a 
ce que les deux conditions suivantes sont satisfaites : 

- le cycle Cs se situe dans la phase d'eff icacite 
decroissante, 

- I'eff icacite au cycle seuil est inferieure a la valeur 
seuil d'eff icacite Es (par exemple Es = 0,1. Eq) : 

Ecs ^ Es 

On peut deja eliminer les points NEC pour lesquels 
I'efficacite est inferieure a Eg. 
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A I'etape 73, on rnodelise la courbe de fluorescence emise 
ajustee dans la decroissance de 1 * ef f icacite, sur la 
portion de cycles Ceep = (Cs - 5) a Cs, grace a I'expression 
(8) dans laquelle on a fait I'hypothese sur la 
5 compensation e telle que s=F'o ♦ 

F„=F'o(l+E'o)'^'=- -POV. -F'o 

Ensuite, le test 74 sur la valeur estimee de E'o et la 

decrementation a I'etape 75 de la valeur du cycle de 
10 basculement Cee? visent a diminuer la valeur recherchee de 
Cee? d'un pas P (qui peut etre egal a 1) et a repeter 

I'etape 73 tant que la valeur de E\ est inferieure a 1. 

Ensuite, lorsque la valeur d'efficacite estimee depasse la 
15 valeur 1 (fleche n en sortie du test 74), on incremente la 
valeur de I'indice Ceep d'un pas h (qui peut etre une 

fraction d'unite) a I'etape 76 et on procede a une etape 
77 de modelisation de la fluorescence Fn du type de celle 

menee a I'etape 73. Tant que I'efficacite estimee E'q est 

20 superieure ou egale a 1 au test 78, on repete les etapes 
76 a 78. Lorsque I'efficacite estimee prend une valeur 
inferieure a 1, on conserve les parametres estimes 

(f^E'o, P\C^Pp) a I'etape de fin 79. 

25 A cette etape, on a obtenu finalement une valeur 

seulement representative de I'effectif initial No dans 
1 ' echantillon d'interet. On peut alors avoir recours a au 
moins un echantillon standard d'effectif initial Nq^^ 
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connu, pour determiner, a I'etape 80, I'effectif initial 
No dans 1 ' echant il Ion d'interet- 



A cet effet, on obtient une valeur de mesure d'un effectif 
5 initial Fost dans un echantillon standard d' effectif 
initial connu Nqs-. On determine ensuite la valeur de 
I'effectif initial dcins 1 ' echantillon d'interet No par 
application d'une meme relation de proper rionnalite entre 
I'effectif initial Nq et sa mesure F'q, dans 1 ' echantillon 
10 standard d'une part, et dans 1 ' echantillon d'interet 
d'autre part. 



En d'autres termes, a I'etape 80 de la figure 9, on peut 
determiner la valeur No de I'effectif initial dans 
15 1 • echantillon d'interet par une simple regie de trois du 

type : Nq = F^^ (n,,^^ / F^^J , impliquant I'effectif initial dans 

le standard Nost et un rapport des fluorescences corrigees, 
compensees et estimees par ajustement du modele de 
fluorescence, pour 1 ' echantillon d'interet et pour 
20 1 ' echantillon standard. 



On comprendra ainsi qu'un seul standard devrait etre 
suffisant pour determiner I'effectif initial de la 
population d'interet dans 1 * echantillon d'interet, selon 
25 un avantage que procure 1' invention. 

Neanmoins, dans une variante qui peut aussi etre prevue si 
necessaire, on obtient les valeurs respectives de mesures 
d'effectifs initiaux F^^, dans plusieurs echantillons 
30 standard d' effectif s initiaux connus Nost- On etablit 
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ensuite une dependance entre les effectifs initiaux Nost 
des echantillons standard et leurs valeurs respectives de 

mesure d'effectif initial F^q^i • Ensuite, apres la 
determination de la valeur de mesure de I'effectif initial 
5 pour 1 * echantillon d'interet F q , on determine I'effectif 
initial de la population d'interet No par interpolation 
sur cette dependance. On comprendra que cette dependance 
peut etre encore typiquement une regression du type 
represente sur la figure 5, mais sur laquelle on reportera 

10 en ordonnees les fluorescences initiales F'q^j des standards 

et de 1' echantillon d' interet F'^ (ou leurs valeurs 
logarithmes respectives), plutot que leur valeur d' indice 
de basculement Ceep- 

15 Des lors que I'on fait appel a un ou plusieurs standards, 
on peut prevoir un ou plusieurs echantillons standard 
d'effectifs initiaux respectifs connus Nostf auxquels on 
applique la succession de reactions d ' amplifications , 
sensiblement dans les memes conditions que pour 

20 1 ' echantillon d'interet. Ensuite, la determination de leur 
valeur de mesure d'effectif initial F^j^^ se fait par 
comparaison des variations parametrees avec les mesures 
experimentales, comme pour 1' echantillon d'interet. 

25 En d'autres termes, on applique bien entendu les memes 
calculs pour les quantites de fluorescence mesurees et 
ajustees sur le ou les standard(s) et pour 1' echantillon 

d'interet. On estime la quantite de fluorescence F^^, avant 

le premier cycle, pour le ou les standard (s) selon le meme 



1 er depot 



33 

precede de determination de pour I'echantillon 

d'interet, decrit ci~avant. 

Un troisieme mode de realisation, correspondant a une 
5 variante du second mode de realisation ci-avant consiste 
globalement a ajuster le modele de I'efficacite En par 
rapport aux mesures experimentales et a injecter ensuite 
ce modele d'efficacite ajuste dans le modele de I'effectif 
courant Nn ou encore dans le modele des mesures Fn. Ce 
10 troisieme mode de realisation se resume comme suit. 

La variation parametree construite a I'etape b) est 
representative du rendement, et, a I'etape c) , on 
determine une variation experimentale du rendement a 
15 partir des mesures experimentales, pour comparer la 
variation parametree a la variation experimentale. 
Ensuite, pour acceder a un parametre representatif de 
I'effectif initial No, a I'etape d) : 

dl) on determine une seconde variation parametree, 
20 representative de I'effectif courant Nn dans 

I'echantillon d'interet, et mettant en jeu au moins le 
parametre representant I'indice de basculement Ceep, 
ainsi qu ' un parametre representatif de I'effectif 
initial Nq, 

25 d2) on applique a cette seconde variation une valeur du 
parametre de I'indice de basculement Ceep determine a 
1 ' etape c) , 

d3) et on ajuste au moins le parametre representatif de 
I'effectif initial No par comparaison directe de la 
30 seconde variation avec les mesures experimentales. 
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On peut avantageusement : 

- a 1 ' etape d2), appliquer une valeur grossiere de 
I'indice de basculement CEEPr comme on a decrit sa 
detection en reference a la figure 4B ci-avant, tandis 
que 

- a I'etape d3), on peut affiner ensuite la valeur de 
1 ' indice con j ointenient avec I'ajustement du parametre 
representatif de I'effectif initial No- 



10 Finalement, on retiendra simplement que le second mode de 
realisation des figures 7 et 8, actuellement prefere, se 
distingue de ce troisieme mode de realisation par le fait 
que 1 ' on ne cherche pas a mener les correlations sur 
1 * ef f icacite, mais on utilise simplement le modele 

15 mathematique de variation de I'efficacite pour modeliser 
et affiner 1' estimation de la fluorescence corrigee et 
compensee . 



Bien entendu, la presente invention ne se limite pas a la 
20 forme de realisation decrite ci-avant a titre d'exemple ; 
elle s'etend a d'autres variantes. 



Ainsi, on comprendra que la presente invention peut 
s' appliquer en outre a la quantification relative, 
25 notamment par PGR. Dans cette application, on prevoit 
d'amplifier, outre la population d'interet, une population 
de reference soit con j ointement dans le meme support, soit 
separement . On releve respect ivement : 

- les mesures exper imentales representatives de 
30 l*effectif de la population d'interet, et 



1 er depot 



35 

- les mesures experimentales representatives de 
I'effectif de la population de reference* 

Le precede peut se poursuivre alors par 1 ' application, a 
la population de reference, des etapes a) , b) et c) , 
tandis que I'etape d) consiste simplement a determiner un 
ratio enLre les effectifs initiaux respectifs de la 
population d'interet et de la population de reference. 

On indique que la quantification relative peut etre 
utilisee pour analyser I'expression d'un gene d'interet au 
cours du developpement d'un organisme. Pour corriger 
notamment les variations de quantite et de qualite entre 
les echantillons preleves sur I'organisrae a differents 
instants, on analyse, en plus du gene cible d'interet, un 
gene de reference connu pour avoir un niveau d' expression 
stable au cours du developpement • 

N 

II s'aqira finalement de comparer les ratios — ^^^^ entre 
les differents echantillons preleves. 

Pour aboutir aux resultats souhaites, deux strategies sont 
possibles . 

Celle de I'art anterieur basee sur la detection du cycle 
Ct se deroule habituellement comme suit. On determine, 
pour chaque echantillon preleve a differents instants tO, 

tl, t2, tn, le ratio —7^-, en utilisant au moins un 

standard ( c' est-a-dire un echantillon dont N^., et A^^ 
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sent connus), ce qui revient a faire deux quantifications 
absolues successives suivies d'un calcul de ratio. 



Une autre strategic particulierement avantageuse dans le 
contexte de 1' invention consiste a determiner, pour chaque 
echantillon preleve a differents instants tO, tl, t2,..., tn 



directement le ratio 



grace a la formule 



Ocihle 



Ocihle 



10 Dans cette seconde realisation qui n* utilise finalement 
que le parametre Fo, en combinant la technique de 
1' invention, aucun echantillon standard n'est necessaire, 
de fagon particulierement avantageuse. 

15 On se refere maintenant a la figure 10 sur laquelle on a 
represents une installation pour la mise en oeuvre du 
precede selon 1' invention. Elle comporte un support SUPP 
comportant ici un puits comprenant 1 ' echantillon d'interet 
ECH et un puits comprenant par exemple un echantillon 

20 standard note St. Le support SUPP est enferme dans une 
enceinte ENC par exemple equipee de moyens de chauffage 
(non representes) pour appliquer au standard et a 
1 ' echantillon d'interet une reaction de PGR. 
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Dans 1' example decrit, on prevoit pref erentiellement 
d'appliquer des mesures de quantites de fluorescence 
emises a chaque cycle par le standard St et par 
1 'echantillon d'interet ECH. A cet effet, on introduit 
5 dans les puits un reactif choisi et on eclaire les 
echantillons a partir d'une lampe (par example halogene- 
tungstene) dans le but de xnesurer des quantites 
respectives de fluorescence emanant de 1 * echantillon 
d'interet et de 1 ' echantillon standard, a chaque cycle de 

10 PGR qui sera applique. En outre, un appareil de detection 
de fluorescence comporte par exemple un objectif 11 pour 
collecter la lumiere issue de la fluorescence, ainsi que 
des moyens de comptage de photons 10, par exemple une 
camera CCD, et/ou des photomultiplicateurs, pour compter 

15 la fluorescence emise a chaque cycle de PGR, par 
I'echantillon d'interet et le standard. Ainsi, la 
fluorescence emise par chaque puits est avantageusement 
focalisee par I'objectif 11, puis detectee 

pref erentiellement par une camera CCD 10 connectee a une 

20 carte d ' acquisition 21, par exemple de type PCMCIA prevue 
dans une unite centrale 20 d'un ordinateur. 



25 



Cet ordinateur est alors relie aux moyens de comptage 10 

precites pour recevoir des signaux representatif s de 
mesures de quantite de fluorescence detectee a chaque 

cycle de PCR et traiter ces signaux pour determiner un 

effectif initial de la population d'interet avant le 

premier cycle, par la mise en oeuvre du procede au sens de 
1 * invention . 



30 
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On indique que 1* unite de traitement comporte 
typiquement : 

- une carte d ' acquisition 21 connectee aux moyens de 
comptage 10, 

5 - une memoire de travail 25 (par exemple de type RAM) 
pour le stockage temporaire et le traitement des 
signaux precites, 

- une memoire permanente 24 pour le stockage du produit 
programme d'ordinateur au sens de 1' invention et des 

10 donnees traitees et pretes a exploitation, par exemple 

pour un diagnostic ulterieur, 

- le cas echeant, un lecteur 22 d ' un support memoire tel 
qu ' un CD-ROM, ou autre, sur lequel peut etre stocke 
initialement le produit. programme d ' ordinateur , 

15 - eventuellement , une interface de communication 26 vers 
un site local ou distant (liaison 28), par exemple pour 
la transmission des donnees traitees pour un diagnostic 
sur le patient realise a distance, 

- une interface graphique 27 reliee typiquement a un 
20 moniteur d'affichage 30, 

- et un processeur 23 pour gerer les interactions entre 
ces differents equipements. 

On indique que 1* ordinateur peut comporter aussi des 
25 organes de saisie tels qu ' un clavier 41 et/ou une souris 
42 relies a I'unite centrale 20. 

On retiendra toutefois que I'' installation au sens de 
1' invention comporte globalement ; 
30 - un support SUPP d ' echantillon au moins pour 
1 ' echantillon d'interet. 
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- un premier appareil ENC pour appliquer, au moins a la 
population d'interet dans 1 ' echantillon d'interet, 
ladite succession de reactions d ' amplification, 

- un second appareil 10 pour relever des mesures 
5 representatives de I'effectif courant de la population 

d ' interet , et 

- des moyens irif oriaa tiques 20 propres a recuperer des 
signaux de mesures du second appareil 10 et a mettre en 
oeuvre tout ou partie des etapes du precede au sens de 

to 1 ' invention . 



A cet effet, un produit programme d'ordinateur peut animer 
ces moyens inf ormatiques . Ce programme peut etre stocke 
dans une memoire de 1' unite de traitement 20, ou sur un 

15 support memoire amovible (CD-ROM ou autre) et destine a 
cooperer avec le lecteur de 1' unite de traitement. Le 
programme d'ordinateur au sens de 1' invention comporte 
aiors des instructions pour la mise en oeuvre de tout ou 
partie des etapes du precede selon 1' invention. Par 

20 exemple, I'algorithme de ce programme peut etre represente 
par un organigramme equivalent au schema de la figure 9. 
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Revendications 

1. Procede mis en oeuvre par des moyens inf ormatiques pour 
quantifier, de fagon absolue et/ou relative, une 
population initiale d'acides nucleiques dans un 
echantillon d'interet soumis a une succession 
d " applications d ' une reaction d ' amplification de la 
population, au cours desquelles on releve des mesures 
experimentales representatives d'un effectif courant de la 
population au moins dans 1 ' echantillon d'interet, 

caracterise en ce qu'il comporte les etapes : 

a) prevoir un modele de rendement (En) de la reaction 
d ' amplification en fonction de la succession 
d ' amplifications , ce modele comportant : 

o une phase sensiblement constante (Eq) pour une 
premiere partie des applications de la reaction, 
et 

o une phase non constante pour une seconde partie 
des applications de la reaction, 
les premiere et seconde parties etant separees par une 
region de basculement du rendement entre les phases 
constante et non constante, incluant un indice 
d ' application (Ceep) correspondant sensiblement au 
basculement, 

b) exprimer, en utilisant le modele de rendement, une 
relation faisant intervenir au moins 1' indice de 
basculement et un parametre representant 1' effectif 
initial dans 1 ' echantillon d'interet (No;Fo), 

c) par comparaison avec les mesures experimentales, 
determiner au moins 1' indice de basculement, et, dans 
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une etape d) ulterieure ou conjointe, en deduire 
I'effectif initial dans I'echantillon d'interet (No). 

2, Precede selon la revendication 1, caracterise en ce 
5 que : 

- a 1' etape b) , on exprime, en utilisant le modele de 
rendement, une variaLion parametree en fonction de 
iadite succession d * amplifications, mettant en jeu au 
moins un parametre representant i'indice de basculement 

10 (Ceep), et 

- a 1' etape c) , par comparaison avec lesdites mesures 
expeirimentales , on determine au moins ce parametre 
representant I'indice de basculement (Ceep) . 

15 3. Precede selon la revendication 2, caracterise en ce que 
ladite variation parametree est representative de 
I'effectif courant (Nn) dans 1 ' echantillon d'interet, 
en ce que ladite variation met en jeu en outre un 
parametre representatif de I'effectif initial dans 

20 1 ' echantillon d'interet (No), 

et en ce que, aux etapes c) et d) , on determine 
sensiblement con j ointement les parametres representatif s 
dudit indice d ' application (Ceep) et de I'effectif initial 
(No) dans 1 ' echantillon d'interet. 

25 

4. Precede selon la revendication 2, caracterise en ce que 
ladite variation parametree est representative du 
rendement, et en ce que, a 1' etape c) , on determine une 
variation exper imentale du rendement a partir desdites 
30 mesures experimentales , pour comparer la variation 
parametree a la variation experimentale . 
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5. Precede selon la revendication 4, caracterise en ce 
que, a I'etape d) : 

dl) on determine une seconde variation parametree, 
representative de I'effectif courant (Nn) dans 
I'echantillon d'interet, et mettant en jeu au moins le 
■ parametre representant ■ ledit indice d 'application 
(Ceep) / ainsi qu ' un parametre representatif de 
I'effectif initial (No) dans 1 ' echantillon d'interet, 

d2) on applique a la seconde variation une valeur du 
parametre de 1 ' indice (Ceef) determine a 1 * etape c) , 

d3) et on ajuste au moins le parametre representatif de 
I'effectif initial (No) par comparaison directe de la 
seconde variation avec les mesures experimentales . 

6, Precede selon la revendication 5, caracterise en ce 
que : 

- a I'etape d2), on applique une valeur grossiere de 
1* indice (Ceep)/ tandis que 

- a I'etape d3), on affine la valeur de 1 ' indice 
con j ointement avec I'ajustement du parametre 
representatif de I'effectif initial (No). 

7- Precede selon la revendication 3 ou la revendication 6, 
caracterise en ce que ladite variation parametree ou, 
respectivement , ladite seconde variation parametree : 

- est representative desdites mesures experimentales, et 

- comporte un parametre cor respondant a une valeur de 
mesure (Fq) representative de I'effectif initial, 

et en ce que 1 ' on determine la valeur de mesure de 
I'effectif initial (Fo) par comparaison desdites 
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variations parametrees (Fn) avec les mesures 
experimentales . 

8. Procede selon 1 ' une des revendications precedentes, 
caracterise en ce qu'on applique une etape prealable de 
traitement des mesures experimentales comprenant une 
soustraction d'un bruit de fond de mesures et 
d ' introduction d'une compensation pour tenir compte d ' une 
mesure non nulle (e) representative de I'effectif initial. 

9. Procede selon 1 ' une des revendications 7 et 8, 
caracterise en ce qu ' on obtient une valeur de mesure d'un 
effectif initial (Fo^^) dans un echaniillon standard 
d'effectif initial connu (Nq^*^) , 

et en ce qu'on determine la valeur de 1' effectif initial 
dans 1 ' echantillon d'interet (No) par application d*une 
meme relation de proportionnali te entre 1* effectif initial 
et sa mesure, dans 1 ' echantillon standard d'une part, et 
dans 1 ' echantillon d'interet d 'autre part. 

10. Procede selon 1 ' une des revendications 7 et 8, 
caracterise en ce qu ' on obtient des valeurs respectives de 
mesures d' effectif s initiaux (Fq^^) dans des echantillons 
standards d'effectifs initiaux connus (Nq^^) , 

et en ce que : 

- on etablit une dependance entre les effectifs initiaux 
(No^^l des echantillons standard et leurs valeurs 
respectives de mesure d'effectif initial (Fq^^), 

- et, apres la determination de la valeur de mesure de 
I'effectif initial pour 1 ' echantillon d'interet (Fc), on 
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determine I'effectif initial de la population d'interet 
(No) par interpolation sur ladite dependance. 

11. Precede selon 1 * une des revendicat ions 9 et 10, 
5 caracterise en ce que 1 * on prevoit un ou une pluralite 

d ' echantillons standard d'effectifs initiaux respectifs 
connus (Nq'*'*^) / auxquels on applique ladite* succession 
d ' amplifications, sensiblement dans les memes conditions 
que pour 1 ' echantillon d'interet, ainsi que la 
10 determination de leur valeur de mesure d'effectif initial 
(Fo^^) par comparaison des variations parametrees (Fn^^) 
avec les mesures experimentales . 

12. Precede selon la revendicat ion 3, caracterise en ce 
15 que 1 ' on prevoit une pluralite d ' echantillons standard 

d'effectifs initiaux respectifs connus, auxquels on 
applique ladite succession d ' amplifications , sensiblement 
dans les memes conditions que pour 1 ' echantillon 
d'interet, ainsi que la determination de leurs indices 
20 respectifs selon les etapes a) , b) et c) , 
et en ce que, a I'etape d) : 

- on etablit une dependance entre les effectifs initiaux 
des echantillons standard (No^^) et leurs indices 
(Ceep^*^) / 

25 - et, apres la determination de I'indice pour 
1 ' echantillon d'interet (Ceep) / on determine I'effectif 
initial de la population d'interet (No) par 
interpolation sur ladite dependance. 

30 13. Precede selon I'une des revendications 1 a 8, 
caracterise en ce que, pour une quantification relative, 
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on prevoit, outre la population d'interet, une population 
de reference soumise a une succession d ' applications de 
reaction d ' amplification, et on releve respectivement : 

- les mesures experimentales representatives de 
I'effectif de la population d'interet, et 

- les mesures experimentales representatives de 
I'effectif de la population de reference, 

le procede se poursuivant par 1 ' application, a la 
population de reference, des etapes a) , b) et c) , tandis 
que 1 ' etape d) consiste a determiner un ratio entre les 
effectifs initiaux respectifs de la population d'interet 
et de la population de reference. 

14. Procede selon I'une des revendications precedentes, 
caracterise en ce que : 

- on exprime les mesures experimentales sous la forme 
d'une variation experimentale du rendement en fonction 
de ladite succession d ' applications, 

- on obtient une variation experimentale du rendement en 
fonction de ladite succession d ' applications, 
comportant : 

* une premiere region sensiblement bruitee pour des 

indices d ' applications faibles, 
* suivie d'une seconde region moins bruitee pour des 
indices d * applications plus eleves. 

15. Procede selon la revendication 14, dans lequel ladite 
phase non constante du rendement est decroissante, 
caracterise en ce que : 

- on estime une valeur grossiere de phase constante du 
rendement (Eq) , et 
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- au moins pour la recherche de 1 * indice (Ceep) dans 
ladite region de basculement , on ignore une partie au 
moins des mesures dans ladite seconde region moins 
bruitee et dont le rendement estime est inferieur a une 
valeur seuil, de preference a une fraction de la phase 
constante (Eo) . 

16. Procede selon 1 ' une des revendications 14 et 15, dans 
lequel ladite phase non constante du rendement est 

decroissante , caracterise en ce que 1 ' on repere ladite 
region de basculement, dans le sens des indices 
d ' applications decroissants a partir de ladite seconde 
region moins bruitee, en detectant un indice grossier 
d ' application d ' amplification dont le rendement depasse 
sensiblement une valeur predeterminee (Eo=l). 

17. Procede selon la revendication 16, caracterise en ce 
qu'on affine 1' estimation de la valeur, le cas echeant 
f ractionnaire, dudit indice d ' application (Ceep) dans la 
region de basculement, dans le sens des indices 
d ' applications croissants a partir de 1' indice grossier, 
en detectant un indice d ' application d ' amplification dont 
le rendement vaut environ ladite valeur predeterminee 
(Eo=l) . 

18. Procede selon I'une des revendications precedentes, 

caracterise en ce que ladite phase non constante du 
rendement est modelisee par une exponentielle decroissante 
comportant un parametre de decroissance (P) , et en ce que 
ledit parametre de decroissance (P) est determine, a 
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li^tape c) avec I'indice dans la region de basculement 
(Ceep) / P^r comparaison avec les mesures experimentales . 

19. Precede selon 1 ' une des revendications precedentes, 
dans lequel la reaction d ' amplification est une reaction 
de PGR (''Polymerase Chain Reaction'') en temps reel. 

20. Precede selon I'une des revendications precedentes, 
dans lequel lesdites mesures sont des quantites de 
fluorescence emises . 

21. Installation pour la mise en oeuvre du precede selon 
I'une des revendications precedentes, caracterisee en ce 
qu'elle comporte : 

- un support d ' echantillon au moins pour 1 ' echantillon 
d ' interet , 

- un premier appareil pour appliquer, au moins a la 
population d'interet dans 1 ' echantillon d'interet, 
ladite succession de reactions d ' amplification, 

- un second appareil pour relever des mesures 
representatives de I'effectif courant de la population 
d'interet, et 

- des moyens inf ormatiques propres a recuperer des 
signaux de mesures du second appareil et a mettre en 
oeuvre tout ou partie des etapes du precede selon I'une 
des revendications precedentes. 

22. Produit programme d * ordinateur , destine a etre stocke 
dans une memoire d'une unite de traitement, ou sur un 
support memoire amovible et destine a cooperer avec un 
lecteur de ladite unite de traitement, caracterise en ce 
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qu'il comporte des instructions pour la mise en oeuvre de 
tout ou partie des etapes du procede selon I'une des 
revendications 1 a 20. 
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